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We solve dynamical equations of motion to determine the conditions under which an over-dense 
region in the early universe will lead to collapse to a black hole, starting from horizon crossing of the 
over-dense region to the point of gravitational instability. Here we focus on the sensitivity to QCD 
and electroweak phase transitions. We then solve rate equations to determine the mass distribution 
of black holes in the present universe. A second order phase transition or rapid crossover would 
have significant consequences only if the index of primordial density fluctuations n > 1.25. However, 
a first order transition would lead to a black hole dominated universe for any realistic value of n 
including n = 1. 



Various symmetries of particle physics that are broken 
in the vacuum are restored at sufficiently high tempera- 
tures These include restoration of chiral symmetry 
and deconfinement in the QCD sector at temperatures 
above about 200 MeV and restoration of the sponta- 
neously broken electroweak symmetry at temperatures 
above about 100 GeV. The impact of cosmological phase 
transitions have generally been found to be inconsequen- 
tial. For example, a first order QCD phase transition 
could in principle have affected nucleosynthesis Q, but 
any reasonable estimate or calculation of the length scales 
for inhomogeneitics in isospin or baryon density yield val- 
ues [l|, H[ that are several orders of magnitude too small to 
have observable consequences Similarly, the baryon 
asymmetry of the universe could in principle have been 
generated during an electroweak phase transition [f| , but 
the consensus is that one must go beyond the standard 
model to generate the observed ratio of baryons to pho- 
tons. In this Letter we shall provide theoretical calcu- 
lations which strongly suggest that if either the QCD 
or electroweak transitions were first order, the early uni- 
verse would have been dominated by black holes, and the 
evolution of the universe would have been very different 
from what it is currently believed to have been. 

Cosmological formation of black holes due to density 
fluctuations has a long history going back to Zel'dovich 
and Novikov Q and Hawking Q. Carr and Hawking 
Q estimated that when the over-density Sp/p exceeded 
a critical value proportional to the square of the sound 
speed v 2 within a particle horizon it is susceptible to 
gravitational collapse to a black hole. Jedamzik |9| sub- 
sequently argued that since during a first order QCD 
phase transition the sound speed vanishes this would be 
a very efficient epoch to form black holes. The physi- 
cal reason is that v 2 = dP/dp, where P is pressure and 
p is energy density, and during the finite time interval 
during which the mixed phase exists there is a change 
in p with no change in P, hence no gradient force to 
prevent collapse. Rather than solving the full equations 
of general relativity for time-evolving fluctuations in an 
expanding universe, which are very computationally in- 
tense, we shall follow an approach used by Cardall and 



Fuller [lfjj. This approach allows us to survey a wide 
class of equations of state and initial conditions. 

We track the evolution of individual over-dense re- 
gions from horizon crossing to the point at which they 
stop expanding, referred to as the turnaround point. We 
model the regions as spherical and homogeneous and 
so are characterized by only two quantities: the epoch 
at which they enter the horizon and their density con- 
trast S = Sp/p at this time. We specify the moment at 
which the regions enter the horizon by the average en- 
ergy density ph of the universe at this time. (Quantities 
with an over-bar are the average quantities for the uni- 
verse. Quantities without the bar refer to those within 
the over-dense region. A subscript h denotes those quan- 
tities evaluated at the time of horizon crossing.) Thus 
p h = (1 + 8)p h . 

The metric for the background is that of a flat FRW 
universe 

ds 2 = -dt 2 + R 2 (t)[dr 2 + r 2 (d6 2 + sin 2 6d4> 2 )] , (1) 

where R(t) is the scale factor of the universe which obeys 
the Fricdmann equation 



dR 
dt 



BttG 



p{t)R 2 (t). 



(2) 



To track the evolution of the over-dense regions we apply 
the metric of a closed FRW universe to the over-dense 



regions 
ds 2 = -dT 2 + S 2 {r) 



dr 2 



1 — nr 2 



r 2 {d9 2 + sin 2 



(3) 

where r is the time coordinate and S(r) is the analog 
of the scale factor for the inner region. The Friedmann 
equation for the region is 



dS_ 

dr 



8ttG 



P (t)S 2 (t)-k 



(4) 



We solve for k by matching the inner and outer metrics 
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at horizon crossing [Io| . 



Sh = Rh 
dS\ ( dR 



dr 



dt 



The evolution of the over-dense region is given by 



fdS_ 
\dr 



8irG 



[p(T)S 2 (r)~j5 h RtS}. 



(5) 
(6) 
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The over-dense region will stop expanding when the 
quantity in parenthesis vanishes. This happens when 



p»S» = phR h 8 



(8) 



where the * denotes the quantity evaluated at this mo- 
ment. 

This approach must be supplemented with the Jeans 
condition to determine whether the region will collapse 
to a black hole. If 



, 5** 7T 
Oh kj 



(9) 



where dh is the particle horizon at crossing and kj is the 
Jeans wavenumber, the region will collapse. The particle 
horizon at crossing is found in the standard way 



d h (t) = R(t) 



dt' 

W) 



(10) 



The relativistic Jeans wavenumber is given by 11 1 



/4ttG(1 + 3v 2 s )w 



(11) 



where w = P + p = Ts is the enthalpy. (The ratio 
of baryon density ub to entropy density s is of order 
10~ 9 and so the baryon chemical potential is neglected 
here.) The critical over-density can now be determined 
numerically if the equation of state is given. 

If the speed of sound is constant the equation of state 
is P = v^p. Then the particle horizon distance can be 
solved analytically with the result that the critical over- 
density necessary for collapse is 



S B = 



8tt 2 



3 (l + v 2 s )(l + 3v 2 s ) 



2\3 



(12) 



For radiation v 2 = 1/3 yielding 5 C = tt 2 /12 w 0.822. 
This is roughly consistent with the numerical results of 
Niemeyer and Jedamzik [l^] who obtained S c = 0.70 ± 
0.02 depending on the exact shape of the initial den- 
sity perturbation when it crossed the horizon. Musco, 
Miller and Rezzolla 13J] and Green, Liddle, Malik and 



Sasaki [14| find 5 C = 0.45 ± 0.02 because they only 
consider the growing mode of the perturbation. For 



this reason the critical deltas obtained using the semi- 
analytical approach as implemented here should be di- 
vided by two. Thus for a purely radiation dominated 
universe 5 C = 7r 2 /24 w 0.411, in satisfactory agreement 
with Musco et al. and Green et al. On the other hand 
this factor of two may just be viewed as a phcnomeno- 
logical factor that better normalizes the results of the 
present approach with the more sophisticated numerical 
calculations. Fortunately our principle conclusions con- 
cerning first order phase transitions are not sensitive to 
this factor of two. 

There are many ways to parameterize an equation of 
state. Perhaps the most physically intuitive is to write 
the entropy density as 



s(T) = ^-T 3 N cS (T) 



(13) 



so that N c g (T) represents an effective number of massless 
bosonic degrees of freedom. Using a physically motivated 
parameterization of iV e ff as a function of temperature one 
can then construct the pressure and energy density from 
thermodynamic identities. For a first or second order 
transition we use the following [lj]. 



N cS (T) 



N 2 - A 2 exp • 
Ni + Ai exp ■ 



T-T c 



r-T e 1 



T>T C 
T<T C 



(14) 



For consistency Ni+Ai < N 2 —A 2 ; equality implies a sec- 
ond order transition and inequality implies a first order 
transition. The iVi is the number of degrees of well below 
T c and N 2 is the number well above. If Aj — > this repre- 
sents a bag model type of equation of state. Although we 
have performed calculations for various strengths of first 
order, second order, and rapid crossover transitions we fo- 
cus here on first order. For QCD we assume three flavors 
of massless quarks plus gluons at high temperature and 
three species of massless pions at low temperatures, plus 
leptons and photons in both phases. We choose T c = 170 
MeV and report here on bag model and softened first 
order with Aj = 0.05T C and Ai = 11.125. For the sec- 
ond order transition we choose A\ = 13 and A 2 = 33.375. 
For electroweak we use a bag model-like equation of state 
with all degrees of freedom above T c = 100 GeV massless 
and a latent heat density varying from 0.5T^ to 1.5T^. 

Figure 1 shows the critical value of over-density for 
QCD as represented by a bag model equation of state. 
(Any strong first order equation of state yields a curve 
that looks very similar.) The remarkable feature is the 
tail. The top side of the tail is determined by the con- 
dition that the over-dense region just enters the mixed 
phase when it stops expanding, and the bottom side by 
the condition that the mixed phase is just ending when 
the over-dense region stops expanding. Any over-dense 
region with initial values within the tail will stop expand- 
ing in the mixed phase and therefore collapse according 
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FIG. 1: Critical value of over-density as a function of energy 
density at horizon crossing for QCD. Within the tail the over- 
dense region is in the mixed phase when it stops expanding 
(dS/dr = 0). A radiation-dominated equation of state would 
have S c = 0.822 independent of energy density. 



to the Jeans criterion. (For the bag model there is even 
an analytic expression for the tail.) A second order or 
rapid crossover transition also have some vestige of this 
tail but it does not extend to infinite energy density or 
approach the x-axis since, although the sound speed is 
small over some range of energy density, it never van- 
ishes. 

The probability that a region which crosses the horizon 
(when the horizon mass is Mh) has an over-density S is 
given by 



P(S,M h ) 



1 



2na{M h ) 



exp 



2a 2 {M h ) 



(15) 



The width <r{Mh) is the COBE-normalized variance at 
the horizon mass Mh taken to be [l~5[ 



a{M h ) = 9.5 x 10" 



M h 



\Q 22 M S . 



(16) 



although it should be noted that Bringmann et al. [lfj 
contend that this is an overestimate. Simple models of in- 
flation predict the power spectrum for the perturbations 
to be a power law ~ k n , where k is the wavenumber as- 
sociated with the perturbation and n is referred to as the 
spectral index. 

To find the present number density of black holes we 
must integrate over the history of the universe. The num- 
ber density at time t is 



N{t) = 



dt' 



R(t) 



1 dV h {t') 



V h {t>) dt' 



1 



v h {f) 



dSP(6,M h (t')). 



(17) 



The last factor is the probability that an over-density 
at horizon crossing will lead to collapse, where A(i') is 
the range of over-densities leading to collapse. Working 
backwards, the next factor is the density factor (one per 
horizon crossing), then the horizon crossing rate, the di- 
lution factor from formation to the observational time i, 
all of which are integrated from to time t. To compute 
the mass distribution dN(M,t)/dM we insert the delta 
function 5 (M co n apse (5, Mh(t'j) — M). What this formula 
leaves out is the possibility that an existing black hole 
will find itself inside another collapsing over-dense region 
at some later time. The observed spectrum of black holes 



dN 



observed 



(M, t) dN(M, t) 



dM 



dM 



exp(-E(M,i)) (18) 



where 
B(M, t) 



dt' 



1 dVh(f) 



V h (t>) dt' 



A(t') 

x e(M coUapse (S,M h (t'))-M) 



dSP(5,M h (t')) 
(19) 



It is this observed distribution that we will subsequently 
refer to. 

Figure 2 shows the resulting black hole mass spec- 
trum at the present time for QCD represented by a bag 
model equation of state. The display begins at 10~ 18 
solar masses since smaller black holes would have evap- 
orated via Hawking radiation by now. The sharp fall-off 
at large masses is related to the horizon mass during the 
epoch of the phase transition, although it does depend 
on the spectral index n. Surprisingly, for n in the range 
from 0.95 to 1.2 the spectrum is represented very well by 
the power-law M~ 4 / 3 . (The deviation from this scaling 
law for larger n is due to contributions from 6 C above 
the critical value for a fixed speed of sound.) All strong 
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FIG. 2: Black hole mass spectrum at the present time due 
to a first order QCD phase transition for various values of 
spectral index n, starting with 1.4 (dashed line) in the upper 
left corner. 
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Table I: Contribution of cosmologically produced black holes to the present value of Q. 
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first order phase transitions have very similar features to 
these. 

The contributions of cosmologically produced black 
holes to the present energy density of the universe, ex- 
pressed in terms of canonical f2, are shown in Table I for 
a variety of equations of state. We originally assumed a 
flat universe and implicitly assumed that black hole pro- 
duction would be a perturbation on the expansion of the 
universe. Therefore any entry in the table with a contri- 
bution to f2 greater than 1.0 would indicate a universe 
dominated by black holes. Since this is apparently not 
the case, a spectral index greater than 1.25 can be ruled 
out no matter which equation of state one looks at. This 
is in close agreement with previous studies [H| . WMAP 
data suggest that n is very close to 1.0, perhaps within 
a few percent [13] ■ However, it must be noted that those 
observations probe comoving wavenumbers which are 10 
to 20 orders of magnitude smaller than those of relevance 
here. Even assuming values of n close to 1 any strong 
first order QCD phase transition would have over-closed 
the universe by many orders of magnitude! (Incidentally 
these results illustrate that the factor of two reduction 
in the S c hardly matters for the tails in Figure 1 since 
the width a(Mh) varies much more rapidly with small 
changes in n.) We have reached the same conclusions 
with a first order electroweak phase transition. For a la- 
tent heat density of the contribution to f2 at present is 
about 10 9 for 0.95 < n < 1.20, and scales with the latent 
heat. The implication is that if either the QCD or the 
electroweak transition were of first order then the evolu- 
tion of the early universe would have been dominated by 
black holes, contrary to prevailing knowledge. How the 
universe would have evolved in such cases is beyond the 
scope of this paper. 

Lattice gauge theory simulations of finite temperature 
QCD indicate that the physical strange quark mass is 
probably too large for QCD to undergo a first order phase 
transition, rather, it undergoes a very rapid crossover. 
Analogous studies in the electroweak sector of the stan- 
dard model indicate a first order phase transition only if 
the Higgs mass is less than around 70 GeV. Since it is 
not, according to lack of observation of the Higgs particle 
at the Tevatron and at LEP, there is no first order phase 
transition in that sector either. It is interesting to pon- 
der the question of how the universe would have evolved 
had either the strange quark or the Higgs particle masses 



been smaller by a factor of 2 or 3. More details concern- 
ing these calculations will be presented elsewhere. 
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